Instability Versus Equilibrium Propagation of Laser Beam in Plasma 
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We obtain, for the first time, an analytic theory of the forward stimulated Brillouin scattering 
instability of a spatially and temporally incoherent laser beam, that controls the transition between 
statistical equilibrium and non-equilibrium (unstable) self-focusing regimes of beam propagation. 
The stability boundary may be used as a comprehensive guide for inertial confinement fusion designs. 
Well into the stable regime, an analytic expression for the angular diffusion coefficient is obtained, 
which provides an essential correction to a geometric optic approximation for beam propagation. 
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Laser-plasma interaction has both fundamental inter- 
est and is critical for future experiments on inertial con- 
finement fusion (ICF) at the National Ignition Facility 
(NIF) 1]. NIF's plasma environment, in the indirect drive 
approach to ICF, has hydrodynamic length and time 
scales of roughly millimeters and 10 ns respectively, while 
the laser beams that traverse the plasma, have a trans- 
verse correlation length, Ic, of a few microns, and coher- 
ence time Tc of roughly a few ps. These microscopic fluc- 
tuations induce corresponding small-scale density fluctu- 
ations and one might naively expect that their effect on 
beam propagation to be diffusive provided self-focusing is 
suppressed by small enough Tc, Tc <IC Ic/cs, with Cs the 
speed of sound. However, we find that there is a collec- 
tive regime of the forward stimulated Brillouin scattering 
(FSBS) instability which couples the beam to trans- 
versely propagating low frequency ion acoustic waves. 
The instability has a finite intensity threshold even for 
very small Tc and can cause strong non-equilibrium beam 
propagation (self- focusing) as a result. 

We present for the first time, an analytic theory of 
the FSBS threshold in the small Tc regime. In the sta- 
ble regime, an analytic expression for the beam angular 
diffusion coefficient, D, is obtained to lowest order in 
Tc, which is compared with simulation. D may be used 
to account for the effect of otherwise unresolved density 
fluctuations on beam propagation in a geometric optic 
approximation. This would then be an alternative to a 
wave propagation code j^], that must resolve the beam's 
correlation lengths and time, and therefore is not a prac- 
tical tool for exploring the large parameter space of ICF 
designs. Knowledge of this FSBS threshold may be used 
as a comprehensive guide for ICF designs. The impor- 
tant fundamental conclusion is, for this FSBS instability 
regime, that even very small Tc may not prevent signif- 
icant self-focusing. It places a previously unknown limit 
in the large parameter space of ICF designs. 

We assume that the beam's spatial and temporal co- 
herence are linked as in the induced spatial incoherence 
[3 method, which gives a stochastic boundary condition 
at z = (z is the beam propagation direction ) for the 



various Fourier transform components H, i?, of the elec- 
tric field spatial-temporal envelope, E, 

E{k,z^O,t)^\E{k)\exp i0kW 

(expi[0kW -0k'(t')]) -f^kk-expl^- li-t'l/T'c). (1) 

The amplitudes, |i?(k)|, are chosen to mimic that of ac- 
tual experiments, as in the idealized "top hat" model of 
NIF optics: 

|£'(k)| = const, k < |ii:(k)| = 0, fc > fc„, (2) 

with 1/lc = k„i ~ fco/(2F), F the optic //#, and the 
average intensity, (|/^ = {^|£'p^ = Iq determines the 
constant. At electron densities, ric, small compared to 
critical, Uc, and for F^ 1, E satisfies (2i] 

fco is ~ the laser wavenumber in vacuum. The relative 
density fiuctuation, p = duc/ne, absent plasma fiow and 
thermal fiuctuations which are ignored here, propagates 
acoustically with speed c^: 

mT'H^+P) ^ (^+2i>|-c?v2) ln(l+p) = cy^L 

(4) 

v is an integral operator whose Fourier transform is vkcs, 
where i/ is the Landau damping coefficient. E is in ther- 
mal units defined so that in equilibrium the standard 
p = exp(— /o) — 1 is recovered. The physical validity of 
Eqs. ®, Q as a model of self- focusing in plasma has 
been discussed before Hi, [13. IfneK 

is taken con- 
stant, there are 3 dimensionless parameters for p ^ I: 
v,h = {kQ/kmY{nclnc)Io/v, and fc = k^CgTc. 

Since Eqn. is linear in E, it may be decomposed, at 
any z, into a finite sum, E = E^^ (x, z, t), where each 
term has a typical wavevector nij : Em.{x,z = 0,t) ^ 
exp(zmj • x). Cross terms E^jE^^, mj ^ m^-, in the 

intensity, vary on the times cale Tc so that their effect on 
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the density response, Eq. 0), is suppressed for Tc <C 1 
(see detailed discussion in Similar consideration 

may be applied to general media with slow nonlinear re- 
sponse, including photorefractive media Then the 
rhs of Eq. Q can be approximated as 



c^V^I = c?V2 V IE, 



j 



= c^V^y" dvF{x,v,z,t). (5) 
F(x,v,z,t) = J , 



n' 

xE^^ (x - r/2, z, t)E*^^, (x + r/2, z, t)e'^-'' / {2nf (6) 

is a variant of the Wigner distribution function which 
satisfies, as follows from Eq. (PJ, 

dF ^ OF if r^, ^, . 
oz ax TT^ J L ^ ' 

exp(-2i[v-v'] -x) -p(2[v- v'],z,i) x 

exp (2i[v - v'] • x) 1 F(x, v', z, t)dv' = 0, (7) 

with boundary value F(x, v, z — 0,t) = Fo(v) = [^^(v)!^. 
Here the unit of x is (l/fco)\/"-c/fT-e and that of z is 
{2/ko)nc/ne. Zero density fluctuation, p = dp/dt — 0, 
is an equilibrium solution of (4), ||SJ| and {Tj), whose lin- 
earization admits solutions of the form, 5p • 
X — Lot)^ for real k and lj, with 
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A =koX/kl-^ = 
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Here and below we assume that the principle branches 
of square and cubic roots are always chosen so that the 
branch cut in the complex plane is on the negative axis 
and values of square root and cubic root are positive for 
positive values of their arguments. The real part of A, 
Ar = Re{\) has a maximum, as a function of uj, close 
to resonance, lo — ±fccs[l + 0{v)]. Below we calculate 
all quantities at resonance uj — ±kcs because analytical 
expressions are much simpler in that case. Xr{k) has a 
maximum, Xmax = kl^Xmax/ko > 0, at /c = kmax, 



krnax/hn - 10^7(3 1^ - 2)22/3c-l + 8 - 21/3^ X 

[31/22(1 +/2)l/2]-\ 

ci = -40 + 225/^ - 27/^, 



C= [Ci+ C2 



,1/3 



C2 - -HIo + 4) V 27 - 60/2 _ 8ij4^ 

Modes with k > kcutoff are stable (A,- < 0), 

with kcutoff = fcm/o(l + ^o)~^^'^/'^^ which defines a 
wavenumber-dependent FSBS threshold. 




FIG. 1: Density source power spectrum, A;*|J(fc,tj)p, with 
/o ~ 4.12, u = 0.15, fc ~ 0.033, and fc^z/fco =i 7.9. The 
dashed lines are at lj = ±fcc. . 



As Iq 0, at fixed k, kgXr —k'^/Io, recovering 
the S{z) behavior of density response function Rq^ in 
If km is set to zero, the coherent forward stimulated 
Brillouin scattering (FSBS) convective gain rate j3| is 
recovered in the paraxial wave approximation. Unlike 
the static response, X{k,uj — 0), which is stable 0| for 
all k for small enou gh In . the resonant response remains 
unstable at small k since as Jo — > 0, Xmax 0.024/o 
and kcutoff ^ kmli/2. 

Since the FSBS instability peaks near u) = ±fccs, 
one expects an acoustic-like peak to appear in the 
intensity fluctuation power spectrum, \I{k,uj)\'^, 
for k less than kcutoff as in the simulation 
(//8, Ue/nc = 0.1) results shown in figure 
1. The fraction of power in this acoustic peak. 



/ \I ik,uj)f duj / J^^ \I {k,uj)f duj 

2/ccs/3<|w|<4fccs/3 / 

increases significantly as /q passes through its threshold 
value for a particular fc, as shown in figure 2. There is 
no discernible difference in shape between \E{k,uj, z)\'^ 
at 2; = 0, where it is oc l/[l -I- (wTc)^] , and at finite z, 
for small Tc. 

If Amax < 1, i-e., Iq < 1, then the FSBS growth 
length, 1/Amax: is large compared to the (vacuum) z 
correlation length, oc fco/fc^, and it is found, for small 
Tc, that a quasi-equilibrium is attained: various low or- 
der statistical moments are roughly constant over the 
simulation range once fc^z/fco ^ 5, as seen in figure 
3. A true equilibrium cannot be attained since (fc^) = 
(|Vi?|2)//o grows due to scattering from density fluctua- 
tions as in figure 4. A dimensionless diffusion coefficient, 
D = (fco/fcm)^(fc^), (proportional to the rate of angu- 
lare diffusion) may be extracted from the data of figure 
4 by fitting a smooth curve to (fc^) for 5 < fc^z/fco < 76, 
and evaluating its slope, extrapolated to z — 0. This 
yields a diffusion coefficient of 4.4E-04. 
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FIG. 2: Fractional power in acoustic peak of the intensity 
fluctuatton spectrum, with parameters as in figure 1, except 
k^z/ko ~ 5.2. Note that the FSBS intensity threshold for 
k/km = 1.5 (2.0) is about 3 (4) 




k^z/k^ 



FIG. 3: A quasi-equilibrium is attained with one point E fluc- 
tuations remaining nearly Gaussian, as evidenced by the small 
change in P5 Q, the fraction of power with intensity at least 
5/0, but strongly modified I — p correlations. Parameters are 
Jo ~ 0.53, f = 0.3 and Tc ~ 0.26. Each curve is normalized 
to its value at z=0. 



D may be compared to the solution of the stochastic 
Schroedinger equation (SSE) with a self-consistent 
random potential 0, p, whose covariance, C^^ (C'' is 
a quadratic functional of F{k)) is evaluated as follows 
[iTJ. Take E as given by Eqn. (|3Jl with p set to since 
it goes to zero with Tc, and use it in Eqn. with 
ln(l + p) ^ p, to evaluate C^. This is consistent only 
if /q < 1, so that the density responce is stable except at 



FIG. 4: For parameters of figure 3, (fc^) = (|V£p)//o, in- 
creases little over the initial equilibration distance of roughly 
5 in these units. The subsequent difi'usion rate is 4.4E-04. 
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FIG. 5: Top hat boundary condition, dashed line, changes 
qualitatively over the propagation distance shown in figure 4: 
solid line at k^z/ko ~ 76. 



small k/k,n- It follows, to leading order in Tc, that the 
SSE prediction for D, for the top hat spectrum. 



DssE = yfJl/&^.^.. 



(10) 



has the value 3.2E-04 for the parameters of Fig. 4. Note 
that DssE is proportional to (p^) and the roughly 20% 
increase of {p^) over its perturbative evaluation (see fig- 
ure 3) used in the SSE accounts for about 1/2 of the 
difference between D and Dsse- 

We find that D depends essentially on the spectral 
form, {\E{k)\'^) ~ F{k) , e.g., for Gaussian F{k) with the 
same value of (fc^), Dcaussam ~ "iD top hat- A numerical 
example of this dependence is found in figures 4 and 5. 
D changes by 40% over 5 < k^z/ko < 76, because F{k) 
changes significantly as seen in figure 5. In this sense, 
for NIF relevant boundary conditions, angular diffusion 
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FIG. 6: Beam angular divergence rate increases rapidly with 
Jo. kl^z/ko ~ 15.7, v = 0.0375, Tc ~ 0.125. In contrast, Eq. 
IIOII predicts a flat curve around 1. 



is an essential correction to the geometrical optics model, 
which (absent refraction) has constant F{k). 

Eqn. (10) implies that d/dz ({k/km)^^ oc 1/km, while 

Amax l/fc^. If the diffusion length is smaller than the 
FSBS growth length, then propagation, which effectively 
increases fc^, will reinforce this ordering. This stability 
condition may be expressed as Z? > Amaxj or qualitatively 
as [13 



(11) 



This is a global condition, as opposed to the wavenumber 
dependent threshold, kcutoffilo)- However, even if Eqn. 
(11) is violated, it is not until kcutoff ~ l-5fcm, so that 
the peak of the density fluctuation spectrum is unstable, 
that FSBS has a strong effect. For these larger Iq values 
a quasi-equilibrium is not attained, and it is more useful 
to consider an integral measure, A{{k'^),z) = {k'^){z) — 
(A:^)(0), of the change in beam angular divergence, rather 
than the differential measure, D. A/Jg is shown in Fig. 
6, normalized to unity at Iq = 0.61. 

Note that we have not observed significant departure 
from Gaussian E fluctuations for Iq < 2 for the param- 
eters of figure 6, which is consistent with the absence 
of self-focusing. Therefore in this regime the effect of 
FSBS is benign, and perhaps useful for NIF design pur- 
poses: correlation lengths decrease, at an accelerated 
pace compared to SSE for Jq ~ 1, with z, while elec- 
tric field fluctuations stay nearly Gaussian. As a result 
[l^ . the intensity threshold for other instabilities (e.g., 
backscatter SBS) increases |2^. If Iq > 4, there are large 
non-Gaussian fluctuations of E, which indicates strong 
self-focusing. 

In conclusion, well above the FSBS threshold we ob- 
serve strong self-focusing effects, while well below thresh- 
old beam propagation is diffusive in angle with essential 
corrections to geometric optics. In an intermediate range 



of intensities the rate of angular diffusion increases with 
propagation. In the weak and intermediate regimes, the 
diffusion results in decreasing correlation lengths which 
could be beneficial for NIF. 
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helpful discussions. 
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